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Abstract 

We investigate a novel quantum random walk (QRW) model, possibly useful in 
quantum algorithm implementation, that achieves a quadratically faster diffu- 
sion rate compared to its classical counterpart. 

We evaluate its asymptotic behavior expressed in the form of a limit proba- 
bility distribution of a double horn shape. Questions of robustness and control 
of that limit distribution are addressed by introducing a quantum optical cavity 
in which a resonant Jaynes-Cummings type of interaction between the quantum 
walk coin system realized in the form of a two-level atom and a laser field is 
taking place. Driving the optical cavity by means of the coin-field interaction 
time and the initial quantum coin state, we determine two types of modification 
of the asymptotic behavior of the QRW. In the first one the limit distribution 
is robustly reproduced up to a scaling, while in the second one the quantum 
features of the walk, exemplified by enhanced diffusion rate, are washed out 
and Gaussian asymptotics prevail. 

Verification of these findings in an experimental setup that involves two 
quantum optical cavities that implement the driven QRW and its quantum to 
classical transition is discussed. 
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1 Introduction 



In a quantum random walk (QRW) (J , a topic of intense research activity in the 
field of Quantum Computing and Information during the last years, the statisti- 
cal correlations between the random coin and walker of a classical random walk 
(CRW), are replaced by quantum correlations. This is achieved by appropriate 
quantization of coin-walker systems and by their dynamic interaction that built 
up entanglement 2 between them in the course of time evolution of the walk. 

Some of the main interests in QRW studies have been the effect of entan- 
glement on various asymptotics, on spreading properties, and on hitting and 
mixing times. From the earlier formulations of QRWs [HI 01 j to recent works on 
general graphs 5 , on the line 6 etc. , it has been shown that some surprising 
features use to distinguish the quantum from classical walks, these included fea- 
tures such as non-Gaussian asymptotics^ J^, quadratic speed up in spreading 
rate for walks on line [Hj, exponentially faster hitting time in hvDercubes|ll)|[TT]. 
exponentially faster penetration time of decision trees U21 CHI j breaking of ma- 
jorization ordering and subsequent increase of entropy on average of position 
probabilities during time evolution of the walk^l] e^c. On the experimental 
side of the investigations on QRWs a number of proposals have been put forward 
for their implementation e.g in ion traps |15| . optical lattices [IBj . or in cavity 

QED [n]. 

In the present work we study a novel model of a QRW the so called model, 
that was introduced in ^|(see alsoj^). This is a discrete time homogeneous 
walk with constant unitary evolution operator V'^, which is the square of some 
V operators (see below). Contrary to all other models of QRW, e.g |5]|15|. 
that evolve by increasing powers of the unitary V, i.e V,V^ ... , the present 
homogeneous model is known to have probability distributions at each step that 
do obey the majorization ordering and have a constant increase in their degree 
of mixing (entropy) ^1] . Also this model is exactly solvable both for its finite 
and long time asymptotic probability distributions, and exhibits a quadratically 
faster spreading rate of its position in comparison to CRW 

The aim of this study is to investigate the effect of changes of the quantum 
coin, identified physically with an atom of two levels modeling the head-tail 
states of a random coin, upon the asymptotic probability distribution func- 
tion (pdf) of the walk. This limit pdf is obtained and has a distinct double 
horn shape. The changes on the coin are induced by letting it to cross an 
optical cavity where it interacts with a single mode EM field. The type of 
interaction implemented in the cavity is the resonant Jaynes-Cummings model 
(JCM)I2niEI|I21|17||2Hl- In fact in order to broad the type of changes made 
on the quantum coin, we will solve for the effects of three other solvable vari- 
ation of the resonant JCM, namely the intensity dependent JCM,23j, the two 
photon JCM[23j and the m— photon JCM[2H1- Since a cavity QED realization 
of QRW has already been proposed [Q, we will say that our cavity driving and 
controlling the walk that precedes that one is the first cavity, while the cavity 
implementing the walk is the second cavity. 

The outline of the paper is as follows: in Chapter 2, the QRW, its moments 
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and its limit pdf arc obtained, in Chapter 3, the effect of the first cavity in 
the form of a completely positive trace preserving map (CPTP), acting on the 
coin density matrix is obtained, in Chapter 4, the conditions for which the 
limit pdf survives or gets destroyed due to effects of the driving cavity are 
obtained and analyzed, and finally in the last Chapter a number of conclusions 
are summarized and some of the future extensions of the work are mentioned. 

2 Quantum Random Walk and its Asymptotics 

Lets us consider a quantum walker system with Hilbert space H^^ = span{ \n)\ n G 
Z} and a quantum coin system with space He = span{\+) , |— )}. Suppose that 
the initial state of the walker is given e.g by the density matrix pw = |0) (0| 
while the initial coin state is a general density matrix pc defined in Hq- The 
evolution operator is taken to be 

V = P+Uo<»E+ + P_Uo<»E_, (1) 

a conditional step operator for the walker state vectors where Uq = e'^"^^ is a 
J— rotation matrix, E± \n) = |n ± 1) are the right/left step unitary operators 
in the walker system, and P+, P_ the projections in the coin space along the 
1+) , |— ) directions. Also important are two other operators, the walker position 
operator L, which satisfies the eigenvalue equation L\n) = n\n), and the phase 
operator $ that acts on the Fourier transformed states H^, = span{\(j)) = 
^ Sn6Z s"™"^ I 4' G [0, 27r)} and admits them as eigenvectors i.e $ |(/>) = 
(f)\<f)}. In terms of this operator it is possible to express the step operators as 
E± = 6=*=** . Then a general density matrix for the walker is written as 

pw= / p(<^,<^') |##' (2) 

Jo Jo 

For initial e.g pw = |0) (0| we have that /9(0, ) = 1. To start investigating the 
dynamics of the walk we suppose we have initially a pure coin density matrix 
PC = \c) (c| , and that the first evolution step involves k applications of V before 
tracing out the coin system. This choice specifics the V'^ quantum random 
walk model. Then the once evolved matrix |(/)) {(p \ is |0) {(p \ £y)=(|</)) (0 |) = 
Trciy'^Pc ® \4>) (0 l^^*')) which is written by means of the Kraus generators 
defined as A±(fc, $;c) = (±|y'=($)|c), as follows 

Eyu (10) (<^' I) = A+ {k, c) 10) {^'\A+{k, c)t + A_ {k, c) |0) {^'\A_ {k, c)t 
= {A+ {k, 0; c)A+ {k,cl>';cr+A_{k, cj>; c)A_ (fc, ; c)* ) |0) (0 | 
= A(fc,0,0';c)|0)(0'|. (3) 

This means that the n-step evolution map operates multiplicative on 10) (0 |,viz. 
£yfc(|0) (01) = A{k, 0, ; c)" 10) (0 I, with A{k, ]c) to be referred to as the 
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characteristic function of the walk. Then the walker density matrix evolves as 



2Tr p27T 



^0 



(4) 



Next is possible to evaluate the quantum moment of the walker position operator 
L, for the n - step evolved density matrix 



Tr{L'el,{pw)) 



1 



2^" 



p(0,<^M"(fc, ;c) 



(S) 



In the last equation P, 



(«) 



;(/9vv)|'7t) is the classical probability for the 



walker to be in the position m after n evolution steps ,and {m'')n are the classical 
statistical moments of the walker position. The asymptotic behavior of these 
moments for large n is 



2tt 



h{2(t>;t)d(t) 



d 



Hence ^ converges weakly to h{(j);k,c) 



(6) 

and 6 as- 



sumes the role of a random variable with probability measure ^^^^"^^ . An alter- 
native expression for the function h useful in our further investigations is given 
below 

h{<j,- k, c) - Tr{[a + y(0)taF(0) + ... + V{<l^)^^^-^) aV{^f-^) pc), (7) 

where a = L^qCtsC/o- 

Proceeding to evaluate the asymptotic pdf of the walk we need to resort to 
the concept of the dual of completely positive trace preserving map[2S|- Con- 
sider the dual map Sy^ : B{Hyy) B{Hw), defined on the set of bounded 
operators acting on the walker Hilbert space Hiy, of some given CPTP map 
Eyk : D{Hw) — > D{Hw), operating on the density matrices pw G D{Hw), 
with operator sum realization Syk (pw) = J2i=± ^ii^, cjpw^iik, c)^. This 
dual map is defined to act on bounded operators X g B{Hw), as e^fc(X) = 
J2i=± ^] cyXAi{k, $; c) . By virtue of the last definition the expectation 
value (quantum moments) of the scaled powers of position operator {^Y , eval- 
uated after n steps of the walker which is now in the state p^,'^ = Sykipw)-, 
become — Tr (^p^\^Y^ , or dually is determined by the equation 

{{^Y)„ = {e*yl{{^Y))o = Tr {pw 4'U(^)^)) • Specializing to the model, 
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Figure 1: The double horn shaped limit probability distribution function for the 
QRW of the V'^ model is given in the plot. Also superimposed to it are given the 
occupation probabilities of the walk as evaluated after n — 36 evolution steps. 
The more spread distribution in the graph is the distribution in (50] rescaled 
appropriately for comparison. 

and taking the case where initially pw — |0) (0| , and pc = |c) (c| , with |c) = 
cosx 1+) + i sinx |— ) , wc get the limit 

lim (e7^ (i-r)) = rn^yp^ = -^^dy. (8) 

\ \ n J / Q Jo 27r 7_i TT^l - y2 

The resulting value of the quantum moment is seen in the last equation to be 
given as a statistical moment of the random variable y'^ with respect to the limit 
pdf determined as follows 

P{y)^—L=, -l<y<l. (9) 
TTVl - 

This pdf, which is the same as that obtained in |29j . determines asymptotically 
the occupation probabilities for the scaled position variable of our QRW. In 
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figure ^ we present its graph which has the shape of a double horn peaked at 
the position y = ±1. This is very much in difference with the Gaussian shape 
of the hmit pdf that occurs in a classical random walk, but shares the double 
horn shape with the limit pdf of another model of QRW[7],ISj, although the two 
distributions differ in their exact functional form. Also in the figure we have 
included the position occupation probabilities of the QRW after a large number 
of steps, in order to show their tendency towards their asymptotic values. 

3 QRW Driven by Quantum Optical Cavity 

We consider now an optical cavity where an interaction between a single elec- 
tromagnetic (EM) mode and a single two-level atomic system takes place on 
resonance. Four different types of this interaction will be considered giving 
rise to four respective solvable models. These are the Jaynes-Cummings model 
(JCM) HjcMi the intensity dependent JCM Hi£i_jcm, the two photon JCM 
H2ph~jCM, and the m— photon JCM Hmph-JCM, with the corresponding Hamil- 
tonians given below: 



HjCM 


= —0-3 -1- wa'a H 


h A((T+a -t- (T~a)), 


(10) 


Hid-. JCM 


= — (T3 + uja^ a H 


h \{a+a^/N + cr^aViV + 1), 


(11) 


H2ph-.]CM 


LO ^ 

— — (T3 + cja'a H 


h A((T+a^ -f cr^a^^), and 


(12) 


Hmph-JCM 


= — 0-3 -1- tja'a H 




(13) 



In the above expressions a\ a are the creation, annihilation operators of 
quanta of the EM field satisfying the canonical commutation relation [a, a^] = 
1, while the (t+,(T^ realize the step operators in the atomic "quantum coin" 
system state space, which together with (73 satisfy the Pauli matrix algebra 
[cr+, (7_] = 20-3, [(T3, a±] = ±(T± . Also A is the field- atom coupling constant and 
u! stands for the atomic energy difference which equals the field frequency in the 
case of resonance, as in our case. 

The dynamic symmetry of the above Hamiltonians can be seen by writ- 
ing them as Hk = C + Vk, k =JCM,ID-ICM, 2ph-JCM, mph-JCM, where 
C — ^(73+ LuN, with N = a^a the number operator, is the common free part 
of the Hamiltonians and represents the total "number of excitations" and the 
Vfc obviously identified with the rest part of the Hamiltonian stands for the in- 
teraction between atom and the field. We can verify the following constants of 
motion: [Hk,C] = 0, [Hk, Vk] = 0, [C, Vk] = 0. 

In the interaction picture the unitary evolution operator reads for each of 
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these models Uk{t) = exp{—itVk)- Explicitly we obtain 

/\, / f\ • sin{ XtVaJa) 

cos(ArvaaT) —la — ^ , . — - , 

u^cMit) = I . I, (14) 

-ta^ — ^ . ^ — ^ cosiXtva^a 

V aaT 



U{t)in-jCM = . ^\ 3i„(,«Lt) (15) 

y —t\/Na' — — ^ cos(AtaTa) y 

and for the m— photon case 

/ C0s(AtV^^^) _iamSBi^^^l \ 

u^..-.c.W^i .^^^^ Jg^ , (16, 

which is specialized to the case of two-photonic transition when m = 2. Assum- 
ing the initial field state is the pure state Pf = |/)(/|, and the atomic coin state 
the general density matrix pc, the state of the atomic system after its crossing 
through the first quantum optical cavity is now described by the reduced density 
matrix given below which is obtain by tracing out the field degree of freedom 
namely, 

euipc) = Trj [u {t)(pc ^ p f )U {t)^) . (17) 
Introducing the operators Pij = in the quantum coin space, then the map 

eu{pc)= E P^iPcPki{f\UlUij\f), (18) 

is a positive and trace preserving transformation of the atomic density matrix 
i.c if pc > 0, then eu{pc) > 0, and also Tr{eu{pc)) = Trpc, as is easily shown. 
For the particular case of the JCM with the field been initially in the vacuum 
state i.e |/) = |0), we obtain the map 

euipc) = SopcSl + SipcSl (19) 
where the so called Kraus generators of that map are 

.o(.) - ( ; ) . - ( s^Ao I ) • 

and satisfy the property S^Sq + SlSi = 1. 

The generalization of this result to the case of all four models when the field 
is in some sharp number state |/) = |r) r = 0, 1, 2, leads to the reduced coin 
density matrix 

suipc) = AipcA\ + A2PcA\ + AspcAI (21) 
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with Kraus generators 



_ , cosiXtrj) \ f 

^ ' cos{Xt0) y ' ^ sm{Xtr]) 



where 77 = VFTT,r + 1, ^/(r + l)(r + 2), ^^^^^ and 6* = V?^, r, r(r - 1), 
are the angles for the respective models viz. the JCM, the ID-JCM, 



(r—m)\ 

2-photon, and the m-photon JCM. The trace preservation of this map requires 
that A\Ai + A\A2 + aIa^ = 1. 

4 Cavity Driven QRW Statistics 

The effect of the first optical cavity will be to transform pc = \c) (c| to eu{pc), 
this yields for the particular choice of |c) = cosx|+) +«sinx|— ), the mixed coin 
state 

suiPc) = ^1 + ^ sin(2x) cos ( Xrjt) cos ( X0t) oi 

+^ [cos (2A77t) cos^ X - cos (2A6'i) sin^ x] 0-3. (23) 
In this case the characteristic function is time dependent and reads 

t) = [- COS {2Xtr]) cos^ x + cos {2Xt0) sin^ x] cos(2(/)) 

+ [sin(2x) cos ( Xtr]) cos ( Xte)] sin(2^) . (24) 

We now proceed to evaluate the limit probability distribution function, and 
to this end we rewrite the last equation as h((f); t) = C{t) cos[2(t!) — A(t)], where 
we have introduce the functions A{t) = — cos {2Xtr]) cos^ x + cos {2Xt9) sin^ X) 
B{t) = sin(2x) cos ( Xtrf) cos ( XtO) from which we define the two functions 
C{t) = y/A{tf + B{tf, and tanA(f) = If we now set y = h{(f);t) then 

for all four inverses h~^, i = 1, 2, 3, 4, we get that h'{h~^{y)) = -2sJC{tY - y^. 
This results into the limit pdf which is i) = 5^ j^T^^^r^^, or finally 

The influence of the first cavity is now expressed by the double dependence of 
the distribution, first on the time spent in the cavity i.e the coin-field interaction 
time, and second on the initial coin state by means of the dependence of C(f), 
via A{f) and B(t), on the angle x of the coin vector. The relationship 

r^^r"^^ (26) 

J-i,/T^ J-\cit)\ Vc{tr - y' 
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between the limit distribution of the walk without the driving cavity (eq.Q), 
and the same one with the driving cavity in presence ('ea. (|25|l '). shows that 
when C {t) is not zero the asymptotics of QRW are robust to the changes caused 
by the modified coin system, up to a scaling. The scaling explicitly refers to 
the changes y{t) y{t)/\C{t)l (-1,1) ^ {-\C{t)\ |C(i))|, in the random 
variable and its interval of values respectively. 

Next, in figure ©J we give the graph of eq.jSBJ at various non zero values 
of C(i). 




-1 -0.5 0.5 1 

y 



Figure 2: The limit pdf of fig.CJ, as it has been affected by the presence of 
the optical cavity. The parameters t and x, are chosen so the function C(t)^, 
modifying the hmit pdf (see text), takes the values C(t)^ = {0.95,0.8,0.7,0.6}. 
The respective plots are with solid line, dotted line, dashed line, and dashed- 
dotted line. Apart from scaling, the limit pdf reappears in the asymptotics of the 
walk. This scaling would lend itself to experimental detection of the robustness 
of the asymptotic pdf of the quantum walk. 

The statistical moment derived by the new distribution are now in general 
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time dependent, e.g the first two of them arc found to be 



lim r ^2cl>;t)^ = 0, 

i^oo n Jq Zn 



(27) 




h{24)-t)h{4)]tf 



d(j) _ C{tY 
2^ " 2 



(28) 



namely the first one is zero so that the walk remains unbiased while the standard 
deviation depends on the time spent by the coin in the first cavity. 

To further probe the behavior of the QRW and especially the influences 
upon its asymptotics of the driving cavity, we first note that for i = 0, that is 
in the absence of the first cavity, the function h depends on the initial state of 
the quantum coin system by means of its angle x i-e /i(0; 0) = — cos(2(/) + 2x). 
However this dependence does not show up in the asymptotic regime since its 
probability distribution given in eq.(|Hl ), appears to have a universal charac- 
ter, and is independent from the initial coin state provided it is of the form 
considered here. 

On the other hand as eas. (|24l25|l indicate in the case that there is a first 
cavity present the dependence of function h on the initial coin state survives in 
the asymptotic regime for the JCM and its variations. Indeed by inspection of 
the functions A{t), B{t) and C{t) as given above we see that they do depend 
on the X angle, and this dependence harbors the possibility of controlling the 
asymptotic statistics of the walk. By choosing the initial coin state to have an- 
gles X = {0, f ,7r, ^}, namely to be |c) = -|+>,-j|-)} respectively, 
we enforce B{t) = 0, Vt. If we further choose the interaction time in the first 
cavity so that A{t) = 0, which implies correspondingly the times t = {-^^^j^, 
(2fc+i)7r ^ (2fc+i)7r ^ {2k+i)T! ^ ^ fc € then wc get also C{t) = for those val- 
ues of t. Then the specific relations among the four coin states and interaction 
times as given above result into only two different pairs of coin density matrices 
and interaction times namely, pc — * = '^tx^^j ^'^'^ Pc = 

t — '^—^^—1 for which C{t) = 0, a fact that holds true for all JC models of our 
study. If each of these two conditions occur we say that a resonance condition 
takes place in the first cavity between the field and the two level atom. So we 
see that the resonance condition implies that the standard deviation becomes 
zero, hence the limit pdf collapses, and more importantly that < >„ ~ n, 
so we loose the quadratic diffusion time speed up, characterizing the quantum 
random walk. In such a case the asymptotic behavior of the standard deviation 
agrees with that of a classical random walk. More precisely what happens is 
that in all the above cases, the exiting coin from the first cavity is in the maxi- 
mally classically mixed state pc = \\-^) {+ \ + \ \—) {—\ = ^1- If this coin system 
is used to feed in the second cavity where the V'^ QRW takes place, then the 
final one-step density matrix for the walker system becomes 



Pw £v^{pw) = -^pw + -tEIpwE^^ + -E^pwE!^. 



(29) 
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A comment is on order. If pw is initially a diagonal matrix then so is finally, 
because of the last equation. Hence we really have a classical one-step transition 
that leads to Gaussian statistics for large n, once we normalize L to ^ . This 
implic's that on resonance the walk becomes fully classical. This analysis makes 
obvious the fact that a judicious choice of the initial coin state permits us to tune 
the interaction time in the first cavity where the JCM or some of its variations is 
implemented, so that we have an absolute control not only over the asymptotic 
behavior of QRW, but on the very quantum nature of its performance as well. 
This conclusion makes the quantum optical experimental investigation of this 
idea worthwhile. 

5 Conclusions 

The QRW of the model V'^ exhibits a quadratically enhanced diffusion rate 
compared to the rate of the classical walk, and an interesting and counter intu- 
itive limit probability distribution. This distribution has been investigated in 
the present work with respect to its behavior under variations of the quantum 
coin state that are tailored in a cavity preceding the black box structure where 
the QRW itself is implemented. As the quantum coin is taken to be a two-level 
atom, its state alterations have been induced by letting it interact with a quan- 
tum mode in a JCM type of interaction on resonance. To gain generality in 
addition to the original JCM, three of its versions namely the intensity depen- 
dent, the two-photon and the m— photon JCM have been used. All four models 
give rise to similar alteration of the asymptotics of the walk, a fact that shows 
the generality of the obtained results concerning asymptotics. 

The two main modifications found of the long time limit probabilities of the 
QRW, are parametrized by the state of the coin going into the cavity, and the 
time spent by the coin in the cavity. In the first modification, the tuning of 
these two parameters to the resonance condition leads to fully classical results, 
namely to Gaussian limit distribution. In the second modification, a comple- 
mentary situation prevails in which the double horn shape limit pdf re-emerges 
up to a scaling. These phenomena arc independent of the particular version of 
the resonant JCM used to realize the driving cavity. Their experimental veri- 
fication seems to be feasible by the present cavity QED experimental settings. 
Additional studies about e.g. the role of interaction time variations, the dcco- 
herence due to spontaneous emission of quantum coins, and the statistics of the 
arrival times of quantum coins in the cavity, are necessary, and will be taken 
up elsewhere. 
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